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We study the diagonalization problem of certain Hofstadter-type models through the alge- 
braic Bethe ansatz equation by the algebraic geometry method. When the spectral variables lie 
on a rational curve, we obtain the complete and explicit solutions for models with the rational 
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geometry properties of the Bethe equation on high genus algebraic curves are investigated in 
cooperation with physical considerations of the Hofstadter model. 
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1 Introduction 



The Hofstadter Hamiltonian is defined by 

H Ho f = p(aU + a" 1 !!" 1 ) + v{fiV + /T^ -1 ), (1) 

where f7, V are unitary operators satisfying the Weyl commutation relation for an absolute value 
1 complex number u, UV = ujVJJ ", and a, /?, it, v are parameters with |a| = |/3| = 1, //, f € R. For 
a primitive A^th root of unity u> (= e 27rv/=:T *), i.e., the phase factor $ = P/A" with P relatively 
prime to N, one may assume the Nth. power identity of U,V: U N = V N = 1. There are several 
physical interpretations for the Hofstadter model, especially in solid state physics. A prominent 
one is to consider it as a tight-binding approximation for electrons bound to atomic sites in a 
two-dimensional crystal and in a strong external magnetic field. The history of the Hofstadter 
model can trace back to the work of Peierls 1 2C ] on Bloch electrons in metals with the presence of a 
constant external magnetic field. By the pioneering works of the 50s and 60s jl], [U], |19L 25], the 
role of magnetic translations was found for the Hamiltonion (|l]) , where the Weyl pair of operators , 
aU and (3V, is a discrete version of magnetic translations in the x, y-directions, with the phase <3? of 
their commutation factor oj representing the magnetic flux (per plaquette), and p, v are the hopping 
amplitudes of the system. Subsequently, a systematic study of this 2D lattice model had begun. 
With the discovery of quantum Hall effects, a large number of interesting and important theoretical 
papers on Hamiltonions of the Hofstadter type appeared in the eighties for the quantum mechanical 



interpretation of the Hall conductivity plateaus (see, e.g., [13] and references therein). In 1976 



Hofstadter [16] found the butterfly figure of the spectral band versus the magnetic flux, in which 
a beautiful fractal picture is exhibited. A detailed study of the model and other Hofstadter-type 
models began thereafter on the explanation on the fractal structure through various mathematical 
approaches, such as the semiclassical approximation, WKB-analysis on difference equation, non- 
commutative geometry and others g, 11, 15, 24]. A pedagogical account of this important 



aspect can be found in a vast literature ( e.g., ^] and references therein). 

On the other hand, motivated by the work of Wiegmann and Zabrodin |23[ on the appearance of 



quantum Uqisl^) symmetry in problems of magnetic translation, Faddeev and Kashaev [12] pursued 
the diagonalization problem on the following type of Hamiltonian by the quantum transfer matrix 
method developed by the Leningrad school in the early eighties (cf, for instant, p2[ ): 

H FK = p(aU + a^U- 1 ) + u((3V + P~ X V~ X ) + p(-yW + ^W' 1 ) , 

where U, V, W are unitary operators with the Weyl commutation relations and the A^th power 
identity property, UV = ujVU, VW = uWV, WU = ujUW; U N = V N = W N = 1. With p = 
in Hfk, one has the Hofstadter Hamiltonion (|l]) with the rational flux. Though the physical 
content of the extended Hamiltonion Hfk with one more added operator W has not yet been 
clarified, this general formulation will provide a nice setting in the study of this kind of Hofstadter- 
type Hamiltonions by the quantum inverse scattering method. With the solution of Bethe ansatz 
equation under a certain postulated degree condition (see (5.26) of p2|), the energy expression of 
d) obtained in |3| 

was reproduced in the p = limit computation. Furthermore, a general frame 



work of calculating spectra of the Hamiltonian through the six- vertex model by the algebraic Bethe 



ansatz method was presented in [12]. In this approach, the Bethe ansatz equation is formulated 
through the Baxter vector^] [||, g], visualized on a "spectral" curve associated to the corresponding 
Hofstadter-type model. In general, the spectral curve is a Riemann surface with a very high genus. 
Thus the Bethe equation here can be viewed as a version of the Baxter's T-Q relation on a high 



It is also called as the "Baxter vacuum state" in other literature. Here we follow the terminology used in 
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genus spectral curve. The method relies on a special local monodromy solution of the Yang-Baxter 
equation for the six- vertex i?-matrix. This solution also appeared in the study of chiral Potts model 
0] . For a finite size L, the trace of the L-monodromy matrix gives rise to the transfer matrix acting 

AT 

on the quantum space ® C , where ./V is the order of the rational flux. The transfer matrix is 
composed of a set of commuting operators, one of which is a Hofstadter-like Hamiltonian. For the 
physical consideration, the size L is kept to be a fixed number. While motivated by the Hofstadter 
butterfly spectral figure, the "thermodynamic" limit will be treated in a manner of the rational 
flux's order iV tending to infinity. This process is a different kind of limiting procedure from many in 
quantum integrable systems, such as the XYZ spin chain or its degenerated forms. For the study of 
models in this paper, we are mainly concerned with the problem of diagonalizing the transfer matrix 
of size L = 3 (see [^], or Sects. 6, 7 of the content). The Bethe ansatz equations of the Hofstadter- 
like Hamiltonians arising in such manner were clearly proposed in 12], and the principle could be 
equally applied to a more general, possibly interesting, models. However, the explicit solutions and 
their qualitative nature have not yet been known, even in the situation when the spectral curve is 
a rational one. It is the object of this paper to show that the algebraic geometry method could 
provide an effective tool for a thorough investigation on the mathematical structure of the Bethe 
equation along the line of the quantum transfer matrix scheme in |l2| . In this work, we obtain the 
complete solutions of the Bethe equation for models with the rational spectral curve for L < 3, 
among which a special kind of Hofstadter-type Hamiltonian Hpx is treated. We present a detailed 
and rigorous mathematical derivation of these solutions of Bethe ansatz equations, including the 
one in [12], and further expanding it to all the other sectors. Both the qualitative and quantitative 
properties of the solutions are studied thoroughly. With the understanding of Bethe solutions 
for these Hofstadter-like models of the rational flux, we apply the thermodynamic process which 
enables us to propose the Bethe equation for a generic flux with the desirable solutions from both 
the mathematical and physical considerations. The analysis we make on the solution of the model 
via the Bethe ansatz method has clearly revealed the algebraic geometry character of the Bethe 
equation. We adopt such an interpretation in approaching the diagonalization problem of certain 
quantum integrable models. Accordingly, the Bethe ansatz method for the Hofstadter model (|l]) 
is formulated along this scheme, and the qualitative nature of Bethe solutions is obtained through 
the approach of algebraic curve theory. 

This paper is organized as follows: In Sect. 2, we outline the concept of Baxter vector, a key 
ingredient for the study of the Bethe (ansatz) equation throughout this work. In order to gain 
the conceptual clarity, we shall present the subject from the algebraic geometry aspect. In Sect. 
3, we first recall some results in [12] relevant to our later discussion, fix notations on the transfer 
matrix and Bethe equation, then derive some general qualitative properties about eigenvalues of the 
transform matrix. In the next four sections, we shall consider the case when the spectral data lie on 
a rational curve, and perform the mathematical derivation of the answer, along with the discussion 
of their physical applications. In Sect. 4, we briefly review the basic procedure of limit reduction to 
a rational spectral curve, originally appeared in [12], and present an explicit form of Baxter vector, 
which will be used in later discussions. We present the complete solutions of the Bethe equations 
of all sectors for L < 3 in Sect. 5. In Sect. 6, we explain the "degeneracy" relations between the 
Bethe solutions and the eigenspaces in the quantum space of the transfer matrix. We also mention 
the relationship between our Bethe solutions and some known results by using the Bethe ansatz 
method in literature, and indicate the difference of our approach of the Bethe equation with the 
usual Bethe-ansatz-type technique through certain non-physical solutions obtained by the latter 
method. In Sect. 7, we will address the "thermodynamic" limit problem in the sense that N tends 
to oo, and discuss the Bethe equation for those Hofstadter-like Hamiltonions appeared previously 
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in the context of Sect. 5, but with a generic q. In Sect. 8, we study the Hofstadter Hamiltonian (|l|) 
by starting from the format in [ O ] . Then we go through an algebraic geometry analysis of the high 
genus spectral curve in which solutions of the Bethe equation are represented, and the connection 
of the spectral curves with elliptic curves are then clarified. We obtain a primary understanding of 
Bethe solutions through the spectral curve, and their relation with the Heisenberg algebra. We end 
with the concluding remarks in Sect. 9 with a discussion of the directions of our further inquiry. 
Convention . In this paper, R, C will denote the field of real, complex numbers respectively, 

n N 

and Z the ring of integers. For positive integers N,n, we will denote <8> C the tensor product of 
n-copies of the complex iV-dimensional vector space C^, and Zjy the quotient ring Z/iVZ. 

2 Algebraic Geometry Preliminary and Baxter vector 

In this paper, we shall denote u a primitive iVth root of unity, and q = lu? with q N = (— 1) N+1 , 

N+l 

in particular for odd N, q = u> 2 . Let Z, X be two operators satisfying the Weyl commutation 
relation with the iVth power identity, 

ZX = ujXZ, Z n = X N = I . 

The algebra generated by Z, X is called the Weyl algebra, in which the element ZX will be denoted 
by Y := ZX. The following relations hold, 

XY = u~ 1 YX, YZ = uj- 1 ZY, y w = (-l)^ 1 . 

The canonical irreducible representation of the Weyl algebra is given by the following expressions 
on the iV-dimensional space C^, 

Z : v 1 ► Z (v ) , Z(v) k = q 2k v k , 

X: X(v), X(v) k = , (2) 

Y : dh Y(v), Y{v) k = q 2k v k „i . 

Here a vector v of C N is represented by a sequence of coordinates, v = (wfc) fcg z> with the iV- 
periodic condition, v k = Vk+N- Hence one can consider the index k as an element of Z^v in what 
follows if no confusion could arise. Denote the standard basis of C^, (k\ the dual basis of 
C N * for k € Tin ■ The /c-th component of a vector v of is given by v k = (k\v). In the Weyl 
algebra, we shall consider only the vector subspace spanned by X, Y, Z and the identity /, in which 
we denote the non-zero operator 

<Pa,toj ■= aY - PX -jZ + 61 : C N — > C N . 

The kernel of the above operator, Ker((/9 Q , j/ 3 j7j( 5), depends only on the ratio of the coefficients, i.e., 
the element [a,P, 7, 5] in the projective 3-space P 3 . The non-triviality of Ker(</? Q , )i g )7) 5) defines a 
hypersurface of P 3 , 

T := {[a, (3, 7, 5] G P 3 | Kei(cp a ^ 5 ) + 0} . 
In fact, one has the defining equation of T as follows: 
Lemma 1 The surface T is defined by the equation, 

T: aF + 6 N = P N + T N , [a,P, 7 ,S] € P 3 . 

Furthermore for [a, (3, 7, 5} € T ', KeT(if a ^ nj s) is an one- dimensional subspace of C N generated by 
a vector v = (v m ) with the ratios, v m : v m -\ = (au> m — (3) : {^u) m — 5). 
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Proof. For v £ C^, it is obvious that the criterion of v in Ker(<£> a a ~ §) is described by the 
above ratio relations of v m s. For a non-zero vector v of C^, the periodic relation v m = i> m +jv for 
a non-zero component u m gives rise to the equation of T. □ 

Let v be a basis of Ker(^ Q! for [a, (3, 7,(5] € .F. Note that there are A^ 2 (projective) lines in 
JF, defined by a N — j3 N = 0, equivalently, 7^ — S N = 0, which are labelled by Pj >fe := {auP — /3 = 
7 — = 0} for j, k £ Zn. Outside these lines, the component of v are all non-zero. For elements 
in P}^, one can set v = \k) at [0, 0,u> k ,l], and v = \j — 1) at [1, cj- 7 , 0, 0]; while for the rest of 
elements, the indices in Z^r with the non-zero components of v form a chain from k increasing to 

i-i- 

For the later purpose, we introduce a family of non-homogenous representations of the surface 
T, depending on the parameter h = [a, b, c, d] G P 3 , 

a = £'a , (3 = xb , 7 = — £'£xc , 5 = — £<i , 

where (x, £, £') € C 3 satisfies the equation, 

Sh : - x V = (-O^O^c* - d N ) . 

For a fixed kP 3 , the operator <p a ,p,f,5 corresponding to becomes 

F(x, £, 0(= F (*> f . /»)) := ^ - x&X + ^xcZ - £dJ , (3) 

and we have 

F(x, £ - 1,0 = £, O - £'xcZ + dJ , £, £' - 1) = F(x, £, £') - ZxcZ - aY . 

For an element p = (cc,£,£') of Sh, we shall denote \p) the basis of Ker(y>£/ aa .j, j<A with 

(0|p) = 1, equivalently, \p) is the vector of defined by 

^ = 1 ^ Hp) _ £W m - x6 



(m — l\p) —^(Z'xcu! m — d) 
We shall call \p) the Baxter vector associated to p G Sh [§, ^, |l2| . Then 

F(x,S,g)]p) = 0, 
F(x,£-U';fc)|p> = |r_p)A_(p), (4) 
= -|r + p)A + (p), 

where Aj- are the following (rational) functions of 5^, 

A-(x,£,£) = d-s^c , A+(x,f,0 = — « > 

and r± are the automorphisms of 5^ defined by t±(x, £, £') = (q x, q -1 ^')- 

3 Transfer Matrix and the Bethe Equation 

As in fll2|| , it is known that the following L-operator^] with the operator- valued entries acting on 
the quantum space C^, 

Lh{x) = ( xcZ ^d ) ' X G ° ' 



3 The conventions in this paper are different from the ones used in where X, Y correspond to our Z, X here. 
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possesses the intertwining property of the Yang-Baxter relation, 

R(x/x')(L h (x) (g 1)(1 (g L h {x')) = (1 (g L h (x'))(L h (x) (g l)R(x/x') , (5) 



aux aux 



where the script letter "aux" will always indicate an operation taking on the auxiliary space C 2 , 
R(x) is the matrix of a 2-tensor of the auxiliary space with the following numerical expression, 



R(x) 



( xuj-x- 1 \ 

cuix-x- 1 ) oj-1 

u) — 1 x — x^ 1 

n —---i 



V 



xuj — x I 



By performing the matrix product on auxiliary spaces and the tensor product of quantum spaces, 
one has the L-operator associated to an element h = (ho, . . . , h,L-i) G (P 3 ) l , 



L-l 



L h( x ) = <g) L h 3 {x) ■= L ho (x) <g> L hl (x) ® . . . ® L/ lt _ 1 (x) , 



(6) 



i=o 



which again satisfies the relation (|5|). The entries of LAx) are operators of the quantum space 
L N 

® C , and its trace defines the transfer matrix, 

TAx) = tr aux (L^(x)) . 

Then the commutation relation holds, 

[T n (x),T n (x')} =0 , x,x'eC. 

The transfer matrix T-Ax) can also be computed by changing L^. to L^. via a gauge transformation 
in the following manner, 



L hj (x) = AjL h {x)A~l l , A L :=A , < j < L - 1 



(7) 



Set 



A; 



1 

1 



and denote the corresponding L^Ax) by (x, £j, Cj+i)- With of (||), we have 

Lhj(x,Cj,Cj+i) = 



Fhj (x, ij - 1, fj+i) -F hj (x, £j - 1, - 1) 
Fhj (x, -F h] (x, (j, - 1) 



and 



where 



T n (x) =tT aux (L n (x,i)) , f := 



L-l 



^(s.O := (g)^(x,^',6+i) 



j=0 



^2,1 (^'0 L h;2,2^ X ^) 



a == Co • 
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The existence of Baxter vectors \pj), pj := (x, £j, £j+x) G Sh 3 for all j with the condition £x, = £q, 
imposes the constraint of elements (po, . . . ,pl~i) on the product of surfaces, Yl^Zo ^hji which form 



curve in Ylj = Q ^hj with the coordinates (x, £o> . . . , £l-i) satisfying the relations, 



l 3= 

f N N _ NuN 

L h ■ «j - I" 1 ) 7n „n — -JW > J-0,...,L-1. (8) 



For p = (po, • • • ,Pl-i) G Cr, the Baxter vector |p) is now defined by 

|p) := \Po) ® . • • ® [p£-i) g 4 C w . 

By the definition of k , the Baxter vector of shares the following relations to entries of 
similar to those for Lr(x,£) in (||), 

^l.l^'^b) = I r -P) A -(P)> ^^(^'^Ip) = l T +P) A +b) > H;2,l( x '0\P) = ' 
where A± , r± are the functions and automorphisms of defined by 

L-l 

A_(x,£ , • • • = II ( d J ~ x ^'+i c i) ' 

A + (*,£o,...,£l-i) = if ^ jdj " , (9) 

r± : (x, • • • , ^ ? _1 £o, • • • , g^L-i) • 

Then follows the important relation of the transfer matrix on the Baxter vector of the curve C^, 

T R (x)\p) = |r_p)A_(p) + \T + p)A + (p) , for p G C g . (10) 

As T^(x) are commuting operators for x G C, a common eigenvector (cp\ is a constant vector of 

with an eigenvalue A(x) G C[x]. Define the function Q(p) = (<p\p) of C?, then it satisfies 
the following Bethe equation, 

A(x)Q(p) = Q(r_(p))A_(p) + Q(r+(p))A + (p) , for p G C K . (11) 

We are going to make a detailed study of the above Bethe equation for the rest of this paper. 
Before that, we first derive certain functional properties of the eigenvalue A(x). 

Lemma 2 With the entries L^,. .(x) of L^{x) in (|6j), the following properties hold. 

(I) Under the interchange of operators, ajY <-> dj , bjX <-» CjZ for all j, we have the symmetries 



among L?..{x)s, L n (x) L? (x) , L? 2 (x) L 



h\2A 



(II) The polynomial L^.. .(x) is an even or odd function with the parity (— l) t+ ^ , and its degree 
is equal to 2[ L+1 2 ^ ] — 1 + 5ij. 

Proof. We apply the gauge transformation @ with Aj = Aj+i for all j. When 

Aj ;= 1 ' 
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one obtains the interchange of entries of Lh{x), hence follows (I). For 




the corresponding L>h(x) is equal to x), which implies the parity of Lr.. (x). The determi- 

nation of the degree of Lr. .(x) can be obtained by certain suitable choices of the values of hjS. 
□ 



From the definition of T^(x),A(x), one can easily obtain the following result. 

Proposition 1 The transform matrix T?(x) is an operator- coefficient even polynomial of x with 
degree 2[4], which is invariant under the substitutions in Lemma || (/). The constant term ofT^(x) 
is given by 

To := 2V(0) = Efto 1 «i ® Y + Ilf=o dj ■ (12) 



L 



Subsequently, the polynomial A(x) in ( |TT[ ) zs an even function of degree < 2[-2-] mf/i A(0) = 
<?' 17^=0 a j + iTj=o d j /° r some l - 

[—] 

By the above proposition, we have T^(x) = Ylj=o 2~2jX 2jf with To given by (12). With a further 
study of the expressions of T%js, one can show that they form a commuting family of operators, 
whose proof we won't present here. Instead, as an illustration of this fact and also for the later use 
of this paper, we list below the explicit form of T2 for L = 2, 3 in T?{x) = Tq + x 2 T2, where the 
commutation relation of T2 and To is easily verified, 

L = 2, T 2 = b oCl X ® Z + c hZ ® X ; 

L = 3, T 2 = b a cia 2 X ® Z ®Y + doh^Y <g> X <g) Z + c aio 2 Z (13) 
+ C0M2Z ® X (g) J + d cib 2 I &Z&X + b Q dic 2 X ® I ® Z . 

The above T 2 for L = 3 can be written as the Hofstadter-type Hamiltonian Hpx of Faddeev 
and Kashaev described in Sect. 1 ( for the exact identification, see J12[|) R. In the equation 
(|TT|) , Q(p) is a rational function of with poles. As the functions of C^, ^j +i x N Cj f — dj and 
(Cf+iO-f ~ x N bj f )(—^j)~ N , are the same, by the description of the Baxter vector, the poles of Q(p) 
are contained in the divisor of defined by 

j=Q 



Hence the understanding of the Bethe solutions of (11) relies heavily on the function theory of Cr, 
the algebraic geometry of the curve will play a key role on the complexity of the problem. We shall 
specify the spectral curves in the later discussion. For our purpose, in the rest of this paper we 
shall only consider the situation when N is an odd integer, and denote the integer [^] by M, 

N = 2M + 1 . 



The operators X, Y, Z, S, T, U in 113] are corresponding to Z, X, Y, U, V, W respectively here in this paper. 



s 



4 The Rational Degenerated Bethe Equation 

For the next four sections, the spectral curve will always be the rational curve under the following 
assumption of degeneration, 

dj = q~ l dj , bj = q~ l Cj for j = 0, . . . , L — 1 . (14) 

To make our presentation self-contained, we shall briefly review in this section the general procedure 
of reducing the Bethe equation on to a polynomial equation, originally appeared in [|TJ], and 
present an explicit form of the Baxter vector which will be convenient for the later use. By replacing 
Cj,dj by 3^,1, we may assume dj = 1 for all j. For the convenience of mathematical discussion, 
also suitable for physical applications, we shall assume that the parameters CjS are all generic. The 
solutions for £j-s in @ are given by 

which possess the structure of a finite abelian group. Therefore is the union of disjoint copies of 
the x- (complex) line indexed by this finite group. Instead working on the curve Cr, the following 
r±-invariant subset of will be sufficient for our discussion of the Bethe equation, 

C := {(x,£o,...,&-i) I £o = ■ ■ ■ = ^l-i = q l , I 6 %n} ■ 

The curve C will be identified with P 1 x Zjy: 

C = P 1 xZ N , (x,q l ,...,q l ) <— » (x,l) . 

The Baxter vectors are now labelled by \x, 1) = ®j =0 \x, l)j, where \x, l)j & C N is defined by the 
relations, 

(k-l\x,l) s - (l-xcrf+n) ' keZN - 

We shall use the bold letter k to denote a multi-index vector k = (ko, . . . , ki-i) with kj £ Z; the 
square-length of k is defined by |k| 2 := J2f=o • The component-expression of the vector \x, I) is 
given by 



j=0i=l 1 xc,(? 



The equation (10) takes the form, 

T(x)\x,l) = \q- l x,l- l)A_(x,Z) + \qx,l- l)A+(x,l) , (16) 
where A± are given by 

L-l L-l 1 _ X 2 C 2 



A^(x,l)=Y[(l-xc 3 q l ), A+(x,0=IIy 



j=o j=o - xc i1~ 

As in we introduce the following two functions on the curve C , 



7=o fc=n i=o fr=o j^ 1 
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by which we define the vectors, 

JV-l N-l 

\x) e m = l*> 2n)/ e (x, 2n)co mn , \x)° m = £ \x, 2n + l)/°(x, 2n + l)u mn . (18) 

n=0 n=0 



By computations, the following relations hold for n G Z 



N, 



rU - 2 " } A ± (x,2n)=A ± (x,0) , ^ r A ± (x, 2n - 1) = A±(s,-1) . 



f^q^x, 2n - 1) ^ v ' 7 IV ' ' ' / e (xg ±1 ,2n-2) 
Then (^) becomes the system of equations, 

T(x)\x)) m = \q~ l x)) m D m (x) + |gx)) m D+(x) , m € Zjy , 

where 

|z))m = (I^Om) I 3 ')™) i D m (x) = ^ ^/"^-(-(x o) J 

In what follows, it is convenient to use the notation of shifted factorial, 

(a; a)o = 1, (a; a) n = (1 — a)(l — aa) ■••(! — aa" -1 ), n 6 Z>o . 
By (|15P (0), we have 

r(x,2n)(k|x,2n) = ^^^1^, 

f(x,2n + l)<k|x,2n + l) = g W a TT \ 1)k ^ n+1 . 

Note that each ratio-term in the above right hand side is defined when the lower index k j + n + 1 
is positive, however its value depends only the class modular N. So we shall use the same notation 
for an arbitrary integer index n by defining the value equal to that of any positive representative 
of the class involved in Z^r, and will keep this convention in what follows. By (|l^), we have 

ntz N T=o ^ "W* 

We proceed the diagonalizing procedure on the matrix D^(x) by a gauge transformation using an 
invertible 2x2 matrix U m (x), 

\x)) m >-> \x) m ■= \x)) m U m (x) , D±(x) ' ^ ^((f^x)" 1 !)^ (x)[/ m (x) . 

We shall choose the matrix U m (x) of the form ^] 



- m 



u(gx) -u(gx) 
u(x) —q m u(x) 



5 The form of the gauge transformation matrix here is slightly different from that of (5.15) in |l2jl b y rearrangement 
of the entries. As we are not able to produce the required formulation through the expression in [p2[ , we wonder that 
there might be a misprint there. 
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Then the diagonalizable criterion of D^(x) for all m is equivalent to the following equation of it(x) 



u(ujx) L -i-} 1 — ax 



/=o 1 ~ c i xq 



Note that the right hand side of the above form is equal to j~7g qv j which is the same as A Q ^ 
The resulting expression of becomes 

D-{x) = g m A_(x,-l) ( J _° x ) , D+(x) = g m A+(x,0) (J . 

With the notations 

l^m = (1^)^71) I 2 ')™.) ) {Qm{ x )i Qm(. x )) = { t P\ x )m > 

one has 

|z)m = k)m9~ m «(^) + \x)° m u{x) , [x)~ = |x)^u(gx) - \x)° m q m u(x) . 
The Bethe equation (111]) now takes the form 



L-l L-l 
± <T m A(x)Q± (x) = J] (1 - xejq-^Qt (xq- 1 ) + R (1 + * Ci )Q± (xg) . (21) 
j=o j=o 

5 Solutions of the Rational Bethe Equation 

In this section, we advance the discussion of the last section to obtain the explicit solutions of the 
Bethe equation (|2l|) for L < 3, among which a special case will be shown in the next section to 



coincide with the one in [12 



Lemma 3 The general solutions of the rational function u(x) for (p0|) are given by 

L-l 

u(x) = R(x N ) J[ ( Cj x; , 
where R(x N ) is a rational function of x N . 

Proof. Note that the ratio of any two solutions of (pc|) is a rational function r(x) with the relation 
r{ujx) = r(x), which is equivalent to r(x) = R(x N ) for a rational function R(x N ) of x N . Therefore 
it suffices to show that YlfZo( c j x 'i u )~m+i ^ s a solution of (p0|), which is easily seen by q = u M+1 . 
□ 



By the expressions of (k|x)^, (k|x)^ in (Jl9|) , in order to have the polynomial functions of Qf^(x) in 
(21), we choose the following gauge function u(x) by setting R(x N ) = n^=o (1 — x N c^) 2 in Lemma 



L-l 

i(x) = Y[(l-x N cf)(x Cj q;q 2 ) M . (22) 
j=0 
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The polynomial Qm( x ) nas the degree at most equal to that of u(x), which is (3M + l)L. By 
Proposition [[], one requires the polynomial solutions Q^(x), A(x) of the Bethe equation pi] ) with 
the constraints, 

deg.Q± (x) < (3M + 1)L 

deg.A(x) < 2[§], A(x) = A(-x), A(0) = q l + 1 for some Z. 

Remark. For another choice of the gauge function u(x), only the function Q^{x) differs by a 
multiple of a certain rational function of x , which makes no effect as far as the Bethe equation is 
concerned. □ 



Lemma 4 Let q be a primitive Nth root of unity, k, I be integers with q k + q l £R. Then q k+l = 1. 

Proof. By the odd property of N, 1 is the only real number among the Nth. roots of unity. We 
may assume that l,q k ,q l are three distinct numbers . The following conditions are equivalent, 

q k + q l G R ^ /-g-'eR ^ q k+l - 1 G q l K . 

By interchanging k and / in the above relations, one concludes that (q k+l — 1) G g fc R fl g'R. By 
g fc ^ ±g', g fc R n q'R consists of only the zero element. Hence q k+l = 1. □ 



For the Bethe equation (|2l|), one needs only to consider the plus-part of the equation because 
of the following result. 

Proposition 2 For m € Zjy, we Zioue Q~ (x) = 0, |x)~ = 0, and 

\x)+ = 2q- m \xY m u{qx) = 2\x)° m u(x) . (23) 

Proof. Let Q~ (x) be a non-zero polynomial solution of ( pip for some A(x) with A(0) = q l + 1, 
and write Qmi x ) = xr Qm*( x ) with Q~*(0) 7^ 0. By comparing the ^-coefficients of (^T]), we have 
-g- m A(0)Q-* (0) = (<T r + q')Q-* (0), hence 

-g- m (l + g I ) = g- r + 9 r . 

By Lemma ||, (/' = q 2m , which implies q r = —q^ m , a contradiction to the odd property of the 
integer N. Therefore the only solution Q~ (x) for the negative-part of is the trivial one. Since 

Qm( x ) is °f t ne form (c/?|a;)~ for an eigenvector (ip\ of T{x) in (g) C^*, and all such vectors (</?| form 

a basis of (g) C^*, hence |rc)^ = for all m. Then follows the relation (p3|). □ 



We are going to derive some general properties of the solutions Q^ n {x) of (21). 

Lemma 5 For a polynomial A(x) with A(0) = q l + 1, the necessary and sufficient condition of A(x) 
for the existence of a non-zero polynomial solution Q^ n {x) of the equation ( |2l|) with the eigenvalue 
A(x) is given by the relation q l = q 2m . In this situation, with Q^ n {x) = x r Q^(x) and Q+*(0) ^ 0, 
one has q r = q ±m . 

Proof. Let be a non-zero solution of (f2~I|) and write Q^(x) = x r Q^.*{x) with Q^(0) ^ 0. 

By the same argument as in Proposition ^, we have 

q- m (l + J)=q- r + q r , 
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hence q l = q 2m and q r = q ±rn by Lemma ^. It remains to show the "sufficient" part of the 
statement. For q l = q 2m , by (O) (p3|), one concludes can not be identically zero. We claim 
that the solution in ( |2l|) has non-trivial solutions. Otherwise, this implies that is 

always the zero vector for all x by the same argument as in Proposition 0, hence a contradiction. 
□ 



Proposition 3 Let m be an integer between and M, and Qm(. x )'Q~N-m( x ) ^ e solutions of ( |2l| ) 
for m, N — m respectively which arise from the evaluation of eigenvectors (<p\ ofTAx) on the Baxter 
vector. Then Q+(x), Q%_ m (x) are elements in x m Ylf=oO- ~ x N c^) C[x}. 

Proof. The divisibility of Q^ n {x), Q~j^_ m (x) by x m follows easily from Lemma ||, so only the factor 
IljQjCl ~~ x N c^) remains to be verified. As Qi~(x) is of the form (<£>\x)f for some vector {(p\ in 

(g) C^*, and by (|l~9|) (p3|), it suffices to show that the following divisibility of polynomials, 



(xc,-; u>)m+i | u(gar) j] 



=0 ( XC i ; W )fcj+n+l 



n'" \ M+l \ I / \ l - / ( XC J' 9 ^ ^ 1 )fc 3 '+n+l 
[xcjuj T ; w) M «(a;) I — ; ; . 

By the form of u(x) in (|2^ ) and the relation quj M+l+ ^ = uj'j +1 for j 6 Z, the above relations are 
easily seen. □ 

For our purpose, the functions Q^ n {x) we shall consider are only those arising from eigenvectors of 
the transfer matrix, hence (x) is in the form of Proposition ||. For the rest of this paper, the 
letter m will always denote an integer between and M, 

< m < M . 



We shall conduct our discussion of the plus-part of the equation (21) for the sectors m, N — m 
simultaneously by introducing the polynomials A m (x),Q m (x) via the relation, 

(A m (s), x m l[(l-x N ^)Q m (x)) = (q~ m A(x), Q+(x)), (q m A(x), Q+_ m {x)) . (24) 
3=0 

Then the relations for both the m and N — m sectors are reduced to the following: 

L-l L-l 
A m (x)Q m (x) = q- m J] (1 - xcjq-^Qmixq- 1 ) + q m J] (1 + xcj)Q m (xq) , (25) 
j=0 j=0 

where Q m (x), A m (x) are polynomials with 
deg.Q m (x) < ML — m 

L 

-2 



deg.A m (x) < 2[^], A m (x) = A m (-x), A(0) = q m + 



The general mathematical problem will be the structure of the solution space of the Bethe equation 
(|25| ) for a given positive integer L. First we are going to derive the detailed answer of the Bethe 
solutions for the simplest case L = 1. 

L = 1. We have A m (x) = q m + q~ m , and deg.Q m (x) < M — m. 
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Theorem 1 The solutions Q m (x)s of (Pq) T _ 1 forms an one- dimensional vector space generated by 
the following polynomial of degree M — m, 

3 m+i—1 „—m—i 

B m {x) := 1 + Y (TT — ^)(xc ) j . (26) 

(Note that the coefficients in the above expression are zero for j > M — m.) 
Proof. Write 

M-m 

Qm{x) = f3j{xc y . 

3=0 

Then (^)l=i is equivalent to the following system of equations of fys, 

( q m + q -m _ q -m-j _ q m+j^. = tfn+j-l _ q-m-3)^ ; j £ Z > , (27) 

where (5k is defined to be zero for k not between and M — m. As the values of q m + q~ m — q~ m ~' J — 
qm+j ; q-m-j _g m +j-i are a rj non-zero, the polynomial Q m (x) is determined by (3q, (or equivalently 
@M-m)i through the recursive relations (|7]). With f3 = 1, this provides the basis element B m (x). 
□ 



Corollary 1 The vector space of all polynomial solutions of (25) L=1 , (without the restriction of 
the degree of Q m (x)), is C[x N ]B m (x). 

Proof. By using the fact 

qm+j-l _ q -m-j _ o j = m -m + 1 (mod N) , 

q m + g -m _ g m+j _ q -m-j = q j = Qj N _ 2m (mod N) , 

and the relation (|27|) for those j with M — m + 1 + ZiV < j < (I + l)iV, (Z € Z>o), any solution 
Sfc>o Pj{cox) k of (P5|) r _i must have (3k = except = 0, ...,M — m (mod iV), hence is an 
element of C[x N ]B m (x) by Theorem □ 



For the Bethe equation (|2q ) with L > 1, by the scaling of the variables, 

x i — ^ A~ a; , Cj i — * Acj for A € C* , 

we may assume that the x 3 -coefficients of polynomials, A m (x),Q m (x), are always homogenous 
functions of cq, . . . ,cl-i with the degree j. As (]25) is invariant under permutations of CjS, the 
coefficients of the polynomials of x involved in (£5) depend only on the elementary symmetric 
functions of CjS, 

Sj = c h--- c ij i for j = 1,. . . ,L . 

i\<...<ij 

We shall denote 

d 

Qm(x) = Y ajx j , d := deg.Q m (x) (< LM - m) , 

3=0 

and define ay to be zero for j not between and d. For the rest of this section, we shall only 
consider the case, L = 2, 3. 
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L = 2. We have only two elementary symmetric functions of Cjs: 

Sl = c + Cl , S 2 = C Cl . 



Lemma 6 Let n be an odd positive integer, A be a n x n-matrix with the complex entries aij 
satisfying the relations 



1 -j 



(-l) i +J+ 1 o n _ J . +1)n _ i+1 , for 1 < i,j < n . 



Then A is a degenerated matrix. 

Proof. Write n = 2h + 1. The determinant of A can be expressed by 



det(A) = s g n ( a ) II = H s S n ( CT ) IT(- 1 



A+a{i)+l 



n— <r(i)+l,n— i+1 



i=l 



i=l 



= sgn(a )(-l) /l+1 ^sgn( f T / )n^'0') J = sgn( f T )(-l) /l+1 det(^) , 

j=l 

where the indices a, a' , run through all permutations of {1, . . . ,n}, and o~$ is the one defined by 
<To(i) = n — i + 1. By sgn(ao) = (— l) h , we have det(A) = 0. □ 



In the equation (|25|) jl=2 ? A m (^) is an even polynomial of degree < 2, and deg.Q m (x) < 2M — m. 
By comparing the coefficients of the highest degree of x in (|25|), we have 

A m ( x ) = ( q m + d + q- m - d - 2 )x 2 s 2 + q m + q~ m . 



For k E Z, we define 



Vk = q h + q k - q m - q m , 



h = (q k 1 -q k )si, 

Uk = ( g fc-2 + g-fc _ g»i+d _ g - m -d- 2 ^ 2 _ 

Then (^) ^=2 is equivalent to the system of linear equations of cejS, 

V m+ ja.j + 8 m+ jO!j-l + Um+jCtj-2 = , j £ Z> . 



(28) 



(29) 



In fact, the non-trivial relations in the above system are those j between 1 and d + 1, hence the 
matrix form of (p9|) is given by 



w m+c( <5m+d u m+d 











V o 











V m+ 2 5m+2 Um+2 

v m+ i 5 m+ i J \ a J 



( OL d \ 
OLd-1 



. 



(30) 



Theorem 2 The equation (J25|) T _ 9 /ios a non-trivial solution Qrn(?^) if and only if deg.Qm^) — 
M — m + m! for < m' < M. For each such m! , the eigenvalue A m (x) in (25) i=2 is egita/ 



to A r 



(x) 



(<1 



m — 1 



+ 



,-m'-2- 



a^s 2 + g m + 



l , and the corresponding solutions of Q m (x) 



form an one- dimensional space generated by a polynomial B m ^ m i(x) of degree M — m + m! with 



B m ,m'(0) = 1. ( Here q^ 



M+l 



) 
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Proof. Denote d = deg.Q m (x). Among those VjS of the entries of the square matrix (fjOl), there is at 



most one zero term, which is given by ujv- 



2m 



0. If Q m (0) = 0, this implies Q m (x) 



„N-2m 



with Q m (0) ^ 0, and each coefficient of the polynomial Q^ n {x) is expressed by a Q^ l (0)-multiple 
of a certain polynomial of cq,c\. By setting c% = 0, x N ~ 2m Q* m {x)Q* m {Q)~' 1 gives rise to a solution 
of (p^) fj=1 , which contradicts the conclusion of Corollary [l]. Therefore Q m (0) ^ 0. If d is less 
than M — m, again by setting c\ = 0, the function Q m (x) gives rise to a solution of (25) i=1 of 
degree < M — m, a contradiction to Theorem |l]. Hence d = M — m + m' for < m! < M. It 
remains to show for each such m' , the solutions Q m (x) form an one-dimensional vector space. As 
any non-trivial solution Q m (x) must have Q m (0) ^ 0, this implies the injectivity of the following 
linear functional of the solution space, 

Qm(x) ^ Q m (0) € C . 

So one needs only to show the existence of a non-trivial solution Q m (x), which is equivalent to the 
degeneracy of the square matrix of size d + 1 in the left hand side of ( |30| ) . Write this square matrix 
in the form 

where C is a (d — 2m') x {2m! + l)-matrix, A, B are the tri-diagonal square matrices of the size 
2m! + 1, d — 2m' respectively. The explicit form of A is given by 



_|_ m / ^M+l+m' 



A 



v M+m' 




>M+m' 







U M+T 











V 







v M+m'-l <W+m'-l UM+m'-l 



VM+2-m' 














\ 



&M+2-m' 
VM+l-m' 



UM+2- 



J 



By (|28|) , we have Vj = VN-j,5j = —5N+i-j,Uj = UN+2-j- This implies that the matrix A satisfies 
the condition of Lemma ^, hence det(^4) = 0. Therefore the square matrix (|3l|) has the zero 
determinant. □ 



Remark. By the above theorem, the following data are in one-one correspondence with integers 
m, m! between and M, 

(m,m) < — ► A m>m /(x) < — >B m>m '(x). 

The characterization of B mm i{x) is given as the (unique) polynomial with &eg.B m ^ m i(x) = M — 
m + m! and B m>m /(0) = 1, whose coefficients ctjS satisfy the equation (^) with 

Vk = q k + q~ k - q m - q~ m , 

5k = (q k ~ l -q~ k )s u 5 (32) 

Uk = (q k ~ 2 + q~ k - q m '~2 - q- m '-i)s 2 . 



L = 3. There are three elementary symmetric functions of CjS, 

si = c + ci + c 2 , s 2 = C Cl + ClC 2 + c 2 c , S3 = c cic 3 . 
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A non-trivial solution Q m (x) of (|25|) has the degree d < 3M — m, and A m (x) is of the form 



A m (x) = \ m x 2 + q m + q- m . 
Note that X m is a homogeneous function of CjS of degree 2. For € Z, we define 



(33) 



q k + q k 



{q k-i 

Jfc-2 



-q )si, 

+ q~ k )s 2 , 



h = (q K 
u k = {q k ~ 3 ~ q~ k )s 3 ■ 
The equation (|25|) L=3 is equivalent to the system of linear equations of ctjS, 

w m +jaj + v m+ jatj-i + (S m+j - \ m )oij-2 + u m+j aj- 3 = , jeZ 



>0 



(34) 



The non-trivial relations of the above system are those for j between 1 and d + 3. 

Lemma 7 Let Q m {x) be a non-trivial polynomial solution of (p5|) for some A m (x). Then the 
degree of Q m {x) is equal to 3M — m with Q m (0) / 0. 

Proof. First we note that for j between m + 1 and 3M, the only possible WjS with zero value are 
given by 

WN-m = W m+ N = . 

Let r be the zero multiplicity of Q m (x) at x = 0. If Q m (0) = 0, by (34) we have r = N — 2m or 
N. The polynomial a~ 1 Q m (x) with c\ = C2 = is a non-trivial solution of (p5|) f =1 with the zero 
multiplicity r. By Corollary [l], r must be equal to N and the degree of Q m (x) is at least N+M — m, 
which contradicts our assumption, d < 3M — m. Therefore Q m (0) ^ 0. By the relation of j = d + 3 
in (p4|), we have q m+d = q—m—d—3^ \i ence the only possible values of d are M — m — 1 , 3M — m. If 
d = M — m — 1, by u>j ^ for m+ 1 < j < M, the function Q m ( x ) with ci = C2 = gives rise to a 
non-trivial solution of (p5|) r=1 with the degree < M — m, a contradiction to Theorem |l|. Therefore 
d = 3M - m. □ 



By the above lemma, the " eigenfunction" Q m (x) for an eigenvalue A m (x) is unique up to a non-zero 
constant, hence there is an one-one correspondence between the eigenvalues and eigenstates of the 
Bethe equation ( |34|) l=3 for a given m. By d = 3M — m, the (d + 3)-th relation in the system (|34| ) 
is redundant, hence the matrix form of ( |34| ) becomes 




, d = 3M 



where A is the N x N matrix, 





1 


u 'n- 


-1 





v' N _ 


2 


6 'n- 


2 


U 'n-2 


w' N _ 


-3 


v' N _ 


-3 


°N-3 



A 



\ 



V o 



w' 



#0 J 



(35) 



(36) 
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with the entries defined by 



Si 



q- '-2 + 



2 + q 2 



)S2, 



2 >3 • 



and 5, C are the following matrices, 



B 



C 



I Sm+i — A m 


UM+l 













\ 


VM 


Sm — A m 


UM 







i 




WM-1 


VM-l 


Sm-i - 













































Wm+3 


A m 
















^m+2 %+2 


~~ A m 




V o 









W m+ i 


"rn+i 


/ 


/ o 


W M +1 


\ 










••• 





WM 




1 


=i m 













™M = WM+1 = q 2 + < 


j 2 — qr m — 


T 


V o ••• 




o J 











Note that the coefficient-matrix of the equation (35) is of the size (d + 2) x [d + 1), while there are 
only d+1 variables ays to be solved. The matrix B is equal to the upper-left (M — m + 1) x (M — m) 
submatrix of the square matrix A — X m I, of which the entries oy, 1 < «, j < N, satisfy the relations 



(— l)* +J OiV-j-|-l,JV-i+l- 



Theorem 3 For < m < M, the condition of the eigenvalue A m (x) = A, 



lX 2 + q m + 



m with a 
= 0, 



non-trivial solution Q m {x) of the equation (pop T _% is determined by the solution of det(A — X m 

For each such A m {x), there exists a unique (up to constants) 



where A is the matrix defined by ( |3q ) . 
non-trivial polynomial solution Q m (x) 
Qm(0) + 0. 



L=3> 



and the degree Q m {x) is equal to 3M — m with 



Proof. By Lemma [^, one needs only to show the existence of a non-trivial solution ays of (35) 
for a X m with det(^4 — A m ) = 0. For m = M, it is obvious as there is no matrix B, and C is zero. 
For m < M, with a given A m , there exists a non-trivial vector in the kernel of A — A m , 



«3M 



— m ^ 



(A - K 



\ CtM-m J 

As the UjS appearing in the matrix B are all non-zero, by the fact that the rank of B is at most 
M — m, one can extend the above ctj (M — m < j < 3M — m) to a solution ays of (^). The result 
then follows. □ 



Remark. By the above theorem, the eigenstate Q m (x) is unique for a given A m (x). It implies 
that for each m, the eigenvalues A m (x) and eigenstates Q m (x) of the Bethe equation (p5|)j,=3 are in 
one-one correspondence. Note that these Q m (x) are obtained under the constraint of A m (x) with 
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the form (|33l ) , a conclusion by the analysis of the transfer matrix in Proposition |lj which we will 
refer as the "physical" criterion while comparing the usual Bethe-ansatz-technique in the discussion 
of the next section. □ 

6 The Degeneracy and Physical Solution Discussion of the Bethe 
Ansatz Relation 

In this section, we first discuss the degeneracy relation of eigenspaces of the transform matrix T{x) 

in <S> C^* with respect to the Bethe solutions we obtained in Sect. 5. As before, A(x) denotes the 
eigenvalues of T(x). By Proposition |], the constant term of T(x) is given by 

T = D + 1, D:=q- L ®Y , 
with the eigenvalue A(0), which is of the form q l + 1. For I E Z^v, we denote 

E^ = the eigensubspace of (g) C^* of the operator D with the eigenvalue q l , 

which is of the dimension N . By Lemma ||, for < m < M, the equation ( |25| ) describes the 
relation of A(x) and Qf{x) through @ when A(0) = q 2m + 1 or q 2 ( N -' m ) + 1. For L = 1, T(x) is 
the constant q~ 1 Y + /, and E' x is the eigenspace of Y for the eigenvalue q l+1 . By the evaluation 
at the Baxter vector |x)^_ m respectively, both the spaces, Bf and Ef~ m , give rise to the 

same functional space generated by x m (l — x N c N )B m (x) with B m (x) in Theorem Q. For L = 2, 3, 
the expression ( |T3| ) of T 2 becomes 

L = 2, T 2 =g- 1 c ci(A:(g)Z + Z® AT) , 

L = 3, T 2 =g- 2 (c ciX(g)Z(g)y + ciC2r(g)X«)Z + coC2^(8)y(8)X) (37) 

+ ^ 1 (c CiZ(g)A"(g)/ + CiC2/®Z(g)X + CoC2X(g)/(g)Z) . 

We shall denote Ol the operator algebra generated by the L-tensors of X, Y, Z, I appeared in the 
corresponding expression of T%. Then Ol commutes with D, hence one obtains a (^-representation 
on E l L for each /. 

L = 2. The 02 is a commutative algebra with the generators X (g> Z,Z ® X, and it contains 

2 N 

the element D(= ZX <S> XZ). The 02-representation on <S> C * has the eigenspace decomposition 
indexed by the (X<8>Z, Z(g> Jf)-eigenvalue (q l , qi), or equivalently, the (D, Z(g>X)-eigenvalue (q l , 
where i,j, I are elements in Z^r with the relation q l = qi +l . In fact, the eigenspace is one dimensional 
with the basis {4>j ,ji I defined by 

(^o,h\--=^ E ^ ok+hkf ~ kk '(k,k'\, 

where (jo,ji) is related to by (a;-' , ) = (<7 1 ,<Z J ). The vectors (0j o jili with a;- 70 " 1 "- 71 = q l , form 

2 

a basis of E 2 . The permutation of tensor factors of <g> C N * induces an automorphism of E 2 , which 
interchanges the vectors {4>j ,ji\ an d (<Pji,jo\- The eigenvalues of T2 are given by CQCi(q 1 ^^ 1 + q-*~ l ) 
for j G Ztv, and the corresponding eigenspace is generated by (^jojj and (^jd, where w- 71 = 
cj j ' 0+Ji = </, with the dimension equal to two for all j except j = (M + 1)/. By (p4|), the index 
which corresponds to (m,m') of Theorem [2] is given by the relation, 

(q l ,q j ) = (q 2m ,q m - m '-^) , (q~ 2m , q- m ~ m ' ~h) . 
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With the evaluation at the corresponding Baxter vector, {(pjojil an d {<f>ji,jo\ gi ye r i se to the same 
eigenstate B m ^ m i[x) in Theorem ^. 
L = 3. We have 

qD = (Z®X®I)(X®I®Z)(I® Z®X) . (38) 

With the identification, 

U = rt®I»J, V = D- 1/2 X®I®Z , (39) 

O3 is generated by U, V which satisfy the Weyl relation UV = q 2 VU and the iVth power identity. 
Hence O3 is the Heisenberg algebra and contains D as a central element. By (|37|), qD~T 2 is the 
following Hofstadter-like Hamiltonion, 

c ci(C/ + U- 1 ) + c c 2 (V + V- 1 ) + Cl c 2 (qD 5 / 2 UV + g^IT 5 / 2 ^ 1 */ -1 ) . (40) 

Note that the above Hamiltonian is a special case of the Faddeev-Kashaev Hamiltonian Hpx with 
W = q~ 1 D~ 5 / 2 V~ 1 U~ 1 , a = (3 = 7 = 1. Our conclusion on the sector m = M is equivalent to that 
in [12] as it will become clearer later on. It is known that there is a unique (up to equivalence) 
non-trivial irreducible representation of O3, denoted by C^, which is of the dimension N. For each 
/, E 3 is equivalent to iV-copies of as 03-modules: E 3 ~ N . For < m < M, we consider 
the space E 3 with q l = g ±2m . The evaluation of E 3 on |x)^ m gives rise to a ^-dimensional kernel 
in E 3 . By Theorem ||, there are N polynomial solutions Q m {x) of degree 3M — m of (|25|) T _% for 
each of N distinct eigenvalues A m (x). The ^-dimensional vector space spanned by those Q m (x)s 
realizes the irreducible representation of the Heisenberg algebra O3. 

Now we discuss the relation between the Bethe equation (^5|) and the usual Bethe ansatz 
formulation in literature. For the physical interest, we are going to focus our attention only on the 
case L = 3. For < m < M, a solution Q m {x) in (^)l=3 must have Q m (0) ^ by Theorem |3[ 
We write 

3M-m -. 

»3M-mQm(x) = J[ (x - -) , Z\ G C* . (41) 



1=1 Z l 



By setting x = z l 1 in (p5|) L _ 3 , we obtain the following relation among z/s, called the Bethe ansatz 
equation in the physical literature as in the case of usual integrable Hamiltonian chains, 

2 3M-m _ 

q m+l TT ^ + Cj_ = yr qzi z n ! < / < 3M _ m . (4 2 ) 

f^m-c 3 J^Zi-qZn 
The following lemma is obvious. 

Lemma 8 For a polynomial Q m (x) of the form (|4l|), the Bethe ansatz relation ( ^2[ ) for the roots of 
Qm(x) is equivalent to the divisibility of q~ m Y\ 2 =0 (l — xcjq'^Q^xq -1 ) + q m Y\ 2 =0 (l + xcj)Q m (xq) 
by Qm{x). In this situation, the quotient polynomial T{x) of the latter pair has the deqree at most 
2 with r(0) = q~ m + q n 



,-,in 



As mentioned in the remark of Theorem |3j, the eigenvalue of the Bethe equation (|25|) ^=3 is de- 
scribed by the eigenstate Q m (x), which is a polynomial with the pre-described zeros satisfying 



(42). However, the quotient polynomial T(x) in Lemma | arising from a solution of ( |42j ) might not 
be an eigenvalue A m (x) with the form ( p3| ) of the Bethe equation (^)l=3, i.e., T(x) might not be 
necessarily an even function. A solution of ([O) with a non-even quotient polynomial T(x) will be 
called a "non-physical" one. For the sector m = M, there is no non-physical Bethe ansatz solution 
of (|4^) by the following lemma. 
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Lemma 9 For m = M , the quotient polynomial T(x) in Lemma || associated to a solution of (E2) 
is always of the form Am(x) in (|33|). 

Proof. By Lemma ||, T(x) = J2]=o1j x ^ an d 7o = Q~ M + Q M ■ It suffices to show the vanishing of 
71 . We write ce^lfQiviix) = x 2M + J2j=o~ 1 Pj x ^ with (3$ 7^ 0. By comparing the x-coefficients of 
dH) L= 3, one has 

-TlPo + 7oA = {-q- M ~ x + g M ) S i/3o + {q- M ~ x + 9 M+1 )/?i , 



which implies 71 = 0. □ 



With a further argument following the proof of the above lemma, one can obtain the relations 
(5.26), (5.27)0 of p^| , where the conclusion applies only on the M-sector. In fact, the comparison 
of the x 2 -coefficients of (|25|) L=3 yields the relation, 

Am — [q + q )S2 + (q - q jsiPiPo + w +9 -9 - q )P2Pq , 

an equivalent expression is the following one, 

2M 

Am = (q~ +q~)s 2 + (q* - q~)si ^ z n + (q? + q~ - q* - q—) s ^ j z i z n . (43) 

n=l i<n 

With the substitution, ^ = qic^ 1 ^ = q2c\ x ,p = qzc^ 1 , the expressions (|42"D m= M; ( |43| ) coincide 



with (5.26) (5.27) in 12]. By Lemma the Bethe ansatz relation ([4^) has been shown to be 
equivalent to the Bethe equation (^)l=3 for the M-sector. However, the parallel statement is not 
true for the sector m = M — 1, in fact, there does exist some "non-physical" Bethe ansatz solutions 
of (|42|). An obvious example is given by the following one. By 3M — m = N, the collection of 
inverse of roots of x N —(3 = 0, {(3 ^ 0), forms a solution of the Bethe ansatz equation ([42j), and 
its associated quotient polynomial T(x) in Lemma |8| is given by 

2 2 3 1 3 1 

T(x) = (q~ + q^) + (q~ — q2)s\x + (q~ + q~)s2X 2 , 



which is not an even polynomial, required by the solutions of the equation (25). This shows that 
the inverse of roots of x N — j3 = for f3 7^ provides a "non-physical" solution of Bethe ansatz 
equation (|42|). By the above example, the constraint ( |33|) on the eigenvalue A m (x) should be taken 
into account in the discussion of Bethe ansatz solutions of ( |4^ ) for the spectrum problem of the 
transfer matrix in an arbitrary sector. Such a consideration will become more crucial when the 
problem of thermodynamic flux limit procedure is involved in the next section. 



7 Rational Bethe Equation for a generic q 

In this section, we are going to study the rational Bethe equation for a generic q. In particular for 
\q\ = 1, it is the infinity flux limit discussion, i.e. N — > 00, of the models appeared in the last two 
sections. In the discussion of this section, q will always mean a generic one. By using the formula 
(H), one has the canonical representation on C°° of the Weyl algebra generated by Z, X with the 
relation 

ZX = q 2 XZ , Y := ZX . 
6 The M in our paper is denoted by P in [Q. 
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With h G (P 3 ) L , one defines the transfer matrix TAx) as in the finite N case, and then discusses 
its spectrum A{x). Under the degenerating assumption (14), the method of the previous three 



sections can equally be applied to reduce the diagonalization problem of the transfer matrix to the 
Bethe solutions of the equation (|25|), It is important to notice that the form of ( |25| ) is valid for 
all finite N while keeping the size L fixed. This enforces us to use the same form of the Bethe 
equation for a generic q as the one in the finite N case, so that the formulation becomes compatible 
with the infinity N limiting process. An eigenvalue of the transfer matrix, the same for the Bethe 
equation, should be the generic analogy of A(x) appeared in the first relation of (p4|), now denoted 
by A m (x) = q m A m (x); however, for the eigenstate of the Bethe equation, we will keep the infinity 
A~ version of Q m (x) in ( p5|) . The Bethe equation for a generic q for m £ Z>o is now described by 



L-l L-l 
A m (x)Q m (x) = J] (1 - xCjq-^Q^xq- 1 ) + q 2m Y[(l + x Cj )Q m (xq) , (44) 
i=0 j=0 

where A m (x) is an even polynomial of degree < 2[4] with A m (0) = q 2m + 1, and Q m (x) is a formal 
(power) series of x, i.e., 

Qm(x) = ^2ajX j ,<E C[[x]} . 

j>o 

As in the discussion of the last two sections, we are going to make a similar analysis on solutions 
of © for L < 3. For L = 1, A m (x) = q 2m + 1, and B m {x) in Theorem [l] defines a formal series, 
which becomes the basis of the solution space of (|44"1) l=i- Note that the description of A m (x) is 
consistent with the spectrum of the operator Y . For L = 2, 3, by the expression (^7j) of T2, we have 
the operator algebra Ol as before, in which To is a central element. 

L = 2. The equation ((44[)_l = 2 is equivalent to the linear systems (p9[) of ays, where the coefficients 
of the equations are defined by (32). By Theorem [2], the eigenvalue A m (x) is given by 



Kn,m'{x) = (q m+m ~> + q" "^)x 2 S2 + q 2m + 1 , m' € Z> . 

For m,m' € Z>o, the solutions Q m (x) of (|44|)l=2 form an one-dimensional space generated by an 

element B m ^ m i(x) S C[[x]] with B mjm /(0) = 1. Note that the algebra O2 is commutative, and the 

2 

spectra of the 02-representation <g> Q°°* give rise to the eigenvalues of T(x), which coincide with 
the above A mjjn i{x) for m,m' 6 Z>o- 

L = 3. We consider the operator qD~T2, which is the Hamiltonian (|40|). For the sector m, its 
eigenvalue is given by q l ~ m \ m , where X m is related to the polynomial A m {x) by 

A m {x) = \ m x 2 + q 2m + 1 . 

With \ m := q~ m X m , the A m s form the spectra of the quadric-diagonal square matrix A ( ptf) f or a 
generic q, which is now of the infinity size as M increases to 00. For each A m , the equation (|4^)^ = 3 
of Q(x) is equivalent to the system (|34"|). It is easy to see that there exists a unique solution Q m (x) 
(up to a constant) by the generic property of q. Note that the same conclusion of Q m (x) equally 
holds for an arbitrary complex number A m , including those not in the spectra of A. This means 
that solutions of the equation (|44|).l=3 alone contain, but not sufficiently determine the eigenvalues 
of the transfer matrix T(x). Those A m s not from the spectra of A correspond to the non-physical 
Bethe ansatz solutions of the finite A^ CELS6, cLS discussed in the last section. 
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8 High Genus Curves and the Hofstadter Model 



We are now going back to the general situation in Sect. 3. Note that the values s of the curve 



in (||) are determined by £q and x N . Denote 



y = x N , v = tf. 

The variables (y, rf) defines the curve 

B K : C n (y)if + (A K (y)-D n (y))r,-B K (y) = 0, (45) 
where the functions At, B^, Cr, D? of y are the following matrix elements, 

\ C K (y) ~D K (y) ) ■ H ^ ycf -< j " 

Note that B? is a double cover of y-line, and is a (Zjv) i+1 -branched cover of B?. The automor- 
phisms t± generate a covering transformation group of over B^. In this section, we shall only 
consider the case 

L = 3 , o = d = , b = c = 1 , 
and we assume the variables hi, h 2 to be generic. The expression of T(x) is given by 

T(x) = x 2 (cia 2 X (&Z ®Y + a x b 2 Z ® Y ® X + b x d 2 Z (g> X <g> J + d ± c 2 X ®I®Z), 

equivalently, x~ 2 D~T(x) is equal to the Hofstadter Hamiltonian ([!]) with U, V given by ( p9| ) and 
[A, v, a, j3 related to hi, h 2 by 

ry OIIIO 0111 

\i = qbic x a 2 d 2 , a = q~ bicj d 2 , v = qa\d\b 2 c 2 , (3 = q a x d\b^ c 2 . 
The curve B^ is defined by (|45|) with 

A n (y) = -y 2 (cf a* + d$<$ ) , B z (y) = y(y 2 c?b» + dfd%) , 
C n (y) = y(y 2 b?c? + af<^) , D K (y) = -y 2 « b$ + b?d$) . 

By factoring out the y-component, we consider only the main irreducible component of Br, denoted 
by 

B : (y 2 b? 4 + a»a»)r? + («W + « - cf - df )yr, - (y 2 c?b» + «) = , (46) 

which is a double-cover of y-line with four branched points, hence it defines an elliptic curve. For 
the curve C^, the variables £0 and £1 are related by = £i N . This implies that consists of N 
irreducible components, each one is isomorphic to the same curve W defined by the equations in 
the variable p = (x, ^0,^,2), 

y\) . c-N _ ~S2 a l + x °1 cN _ ~SQ a 2 + x °2 

VV ■ ^0 - NpN r N jN ' ^2 - N tN r N jJV • 
X ?2 L l u l x SO L 2 a 2 

It is easy to see that W is a A r3 -fold (branched) cover of the elliptic curve (|46|), and the genus of 
W is 6./V 3 — QN 2 + 1. We shall label the irreducible components of by s £ Z^r and denote the 
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elements of by (p, s) with p £ W, whose .^-coordinate is given by £o£i = Q 2s ■ The relation 
(10) now becomes 



T(x)\p,s) = \r-(p),8- l)^- i A„(p) + |r+(p), S - iy s A+(p) , 

where r_|- are transformations of of W defined by (^), with the coordinates only involving (x, £0,^2)1 
and A-t are the following functions on W, 



A-O) = -xg l {x£ 2 ci - di)(x£oc 2 - d 2 ) , 

d\— x 2 b\C\){a 2 d 2 — x 2 b2< 
(&ai-xb 1 )(£ a 2 -xb 2 ) 



A — rvt (aid 1 -x 2 b 1 c 1 )(a2d2-x 2 b2C 2 ) _ r c c \ _ yM 



By averaging the vectors |p, s) over an element p of W, one defines the following Baxter vector on 
W 

1 Ar_1 2 

|p) := a? 12 \p> s )<i s2 ■ 

The action of the transform matrix T{x) on \p, s) can be descended to the Baxter vector of W as 
follows: 

x- 2 T{x)\p) = |r_(p))A_(p) + |r+(p))A + (p) , (47) 
where A± are the functions of W, 

(xg2 CI - rfi ) (xgp C2 - (fe ) 



A_(x,e ,&)= ^-°^g 

A (t- £ £ — g2(ai^l-a: 2 feici)(a2rf2-^ 2 t'2C 
^+\ X ^0,C,2) - x{i 2 a 1 -xb 1 )(^a2-xb 2 ) 

By the following component expression of the Baxter vector on C^, 



fj[ ^xci^-dx &{£oxc 2 u>3 - d 2 ) 

the Baxter vector on W is given by 

= ES^ g - 2 ^_ 2fcofcl _ fc!+fel ^ $o(-g 2 a 1 ^ + x6 1 ) |J (48 ) 



iV 



* g 2 xci^* - di 4^ £ 2 (£o^c 2 uP - d 2 ) ' 



Note that the kjs in the above formula are integers modular N. Each product-term of the right 
hand side means the one for a positive integer kj representing its class in Z^r- For an eigenvector 

3 

(<p\ in (8 C^* of the operator x _2 T(x) with the eigenvalue A, by (|47|), the function Q(p) := (<p\p) 
of W satisfies the following Bethe equation, 

\Q(p) = Q( r _(p))A_(p) + Q(r+(p))A + (p) , A € C . (49) 

The above equation possesses a Z 2 -symmetry with respect to the following involution of W, 

<r : W ► W , p = (x, 6) ff(p) = (-x, -Co, -6) • 

In fact, the commutativity of a and t±, and the u-invariant property of A±(p) are easily seen. 
Then by (H8|), Baxter vector \p) is invariant under a, i.e., |p) = |<r(p)) for p G W, which implies 
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that Q(p) is a cr-invariant function. Furthermore the rational function Q(p) has the poles contained 
in the following divisor, 

N-i 

C^" 1 II &x Cl J - diKZoxczJ -d 2 ) = 0, 
i=i 

in particular, it is regular at x = 0, oo. Note that the finite values of Q(p) at x = 0, oo are consistent 
with the asymptotic values of A± at x = 0, oo in the equation (|49|), 

A±(x,^o,6)= ±x-^o + O(l) asx^O, 
A±(x,£o,£2) = ±x&2CiC2 + O(l) as X — > OO . 

With the x-coordinate, W is a 2iV 2 -cover over the x-line, unramified at points over x = 0, 00 whose 
(x, £o, ^-coordinates are given by 

ii , = ±{0,qJ , q \ — — ) , oo fi , = ± {00, q J— -,q W- , 1,1 EZ N . 50) 

Consider the -D-eigenspace decomposition of ® q n * = 0^ g ^ E 3 in Sect. 6. The evaluation on 
the Baxter vector of W gives rise to the following linear transformation, 

ei : E 3 — ► {rational functions of W}, f 1— ► e;(v)(p) := (i>|p) , for Z G Z^r • 



Theorem 4 For Z € Zjy ; i/ie linear map e\ is injective, hence it induces an identification of E 3 
with a N 2 - dimensional functional space ofW. 

Proof. Define the following vectors in (g) C^*, 
where G Zjy. We have 

(i) koMM \Z g> X ® 7 = w *o+ fc i {tpk MM I . 
(tfkoMte \X ® I ® Z = u k2 (ipko-iMte I > 
(^fc ,fc 1 ,fc 2 |-Z'® Z®X = (^ fc0)fcl+ i )fe2 _i| . 

By (|38D, one has 

(^fcoA.tel^ = ^ 1+2(fc0+fcl+fc2) (^ -i,fe 1+ i,A: 2 -il , (^o,h,n\D = q- Wjo (d> jo ,n,n\ ■ 

For a given Z G Z^r, let jo be the element in Zjv defined by q l = g _1+2j0 . Then the vectors (^o,ji,i2l 
with ji , j 2 G Ztv form a basis of E 3 . By ( |48l) , we have 

(^ 0)fel , fe2 b)=iv- 1 ^ e g - fcq +^- 2 *°+ i ) fc> fi * o( ; 6oi f ft rfr^t xh ^ > 
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hence 



"jo 



-3k 2 +k(-2j +6ji+2j 2 +l)-k' 2 +4k'k+(-2j 1 +l)k' 



k' go(-g2ai^'+xfei) y-rj2 -k -gpc^t^+zfe 



i=1 ^xcx^-dx '■'■3=1 



Set p = 0^, defined in (|50|). By the relations of their (£0, ^-coordinates, £o£2<2i^i 1 = , 
?o^2~ la 2^2 1 = we obtain 

jlh \0f.,) = AT-ig-^+^+iaCl+i-i') ^ (? -2fc 2 +fc(-2jo+6ii-i+i')^'(-2ji+2+i+i'+4fc) 

= gJ 0(l- J l)-2 J1 + 4+i2+i2 2 + - Jl(3 ' + ' , l +2j2( '"' ,) + 2 ^- ( ' + ' ,+ ^ 

hence 

^-l+^V+^l-^-- h-i 2 /A lrd= \ ('+' , +2)(-i+3/+2)-4j (»+/) i(-3ji+2j 2 )-i'(j 1+ 2j 2 ) 

As the correspondence, (j'1,,7'2) l— * ( — 3ji + 2j2,— ji — 2j2), defines an automorphism of 7? N , the 



relation (51) gives rise to an isomorphism between and the space of Baxter vectors 0^~-,s (or 



equivalently O^s). This implies the injectivity of □ 

By the discussion in Sect. 6, E3 is equivalent to N copies of the standard representation as the 
Heisenberg algebra 03-modules. Hence by Theorem |], there exists an C?3-module structure on 
£j(E 3 ), inherited from the representation space E 3 . The mathematical structure of the functional 
space q(E 3 ) by incorporating the divisor theory of Riemann surfaces into the Heisenberg algebra 
representation remains an algebraic geometry problem for further study. 

9 Conclusions and Perspectives 



We follow the framework in [|12[ by the quantum integrable method to study the diagonalization 
problem of some Hofstadter-like models. Through the Baxter vector of the spectral curve, the 
study of diagonalizing a Hamiltonian with the rational magnetic flux is reduced to the problem of 
a certain " Strum-Liouville-like" difference equation on the curve, called the Bethe equation of the 
associated model. The spectral curve is in general with a large genus, and the relations among zeros 
and poles for a solution of the Bethe equation yields a system of algebraic equations. Such systems 
of relations among zeros of the Bethe polynomial solution on a rational spectral curve are usually 
referred to the Bethe ansatz equation in literature. For certain models of physical interest, e.g., the 
Hofstadter Hamiltonian ([]]), the study of the high genus spectral curve is a necessary step in solving 
the spectrum problem through the algebraic Bethe-ansatz-technique. In this paper, we clarify some 



finer mathematical manipulations in [12], then go through all the delicate points one must consider 
in order to obtain the explicit Bethe solutions. A careful analysis of the mathematical nature of the 
the Bethe equation reveals the vital role of algebraic geometry in a thorough understanding of the 
Bethe ansatz method, so is the same with the need to obtain the physical answer of the associated 
model. For this reason, we have examined in this paper the Bethe ansatz equation in the context 
of algebraic geometry, even in the degenerated rational spectral curve case in order for gaining 
the mathematical insight of the Bethe equation. We further extend the approach to some more 
general situations. Above all, we have endeavoured to present a clear and self-contained account 
of the theory, and hope to have elucidated the mathematical structure of the Bethe-ansatz-style 
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method in the physical literature. The main content of this paper is in the discussions after Sect. 
5, where the detailed mathematical derivation and analysis comparable to physical considerations 
are presented. The topics between Sect. 4 and Sect. 7 are devoted to the degenerated case, where 
the Bethe equations are related to models with the rational spectral curve. With an explicit gauge 
choice, we have conducted the mathematical investigation of the Bethe equation and obtained the 
complete Bethe solutions for all sectors in Sect. 5. With these mathematical results, we are able to 
further advance the study of the relevant physical problems, namely, the "degeneracy" of eigenstates 
of the transfer matrix to the Bethe solutions and the thermodynamic limit discussion in Sects. 6 
and 7. Furthermore in Sect. 6, the explicit calculation we have performed for the Bethe solutions, 
when specializing on one particular sector, provides the results parallel to those using the usual 
Bethe ansatz method in Meanwhile, the finding of some non-physical Bethe ansatz solutions 
in other sectors supports the justification of our approach of the problem. The method we employ 
here can apply equally well to any number of size L. However, to keep things simple, we restrict 
our attention in this paper only to the case L < 3 in the discussion of Bethe solutions; the analysis 
on L = 1, 2 are mainly for the mathematical purpose to pave the way of discussing the models with 
a higher size L. The Hofstadter-like model is related to L = 3, of which the Bethe ansatz equation 
has been mathematically discussed in details here. For L = 4, it is expected that the problem 
should be closely related to the discrete quantum pendulum by the works of (9|, [L8|]. The results 
we have obtained in this paper through the Bethe equation approach strongly indicate a promising 
direction to the spectra problem of other models, e.g., the discrete quantum pendulum. For the 
thermodynamic limit discussion, the analysis in Sect. 7 on the special Hofstadter-like Hamiltonian 
(|40|) for a generic q shows that the Bethe relation on the spectrum and eigenstates we proposed are 
in accordance with the diophantine approximation process of an irrational flux. The comparison 
of our Bethe ansatz method with the C*-algebra approach of semiclassical analysis employed in || 
for the multifractal spectrum structure is a fascinating problem. We plan to address the question 
of a such program elsewhere. 

For the original Hofstadter model (|l]), the Bethe equation is formulated as a "difference" equa- 
tion of functions on a high genus spectral curve. In Sect. 8, we have made a primary investigation 
on its solutions . As the spectral curve is a Galois abelian cover of an elliptic curve with the 
covering group determined by the order N of the rational flux, it would be essential to have a de- 
tailed algebraic geometry study of such a high genus curve in cooperation with the Bethe solutions, 
so that the base elliptic curve and the classical elliptic functions could be engaged in the theory. 
With our finding in the rational spectral curve case as a guidance for some appropriate direction 
of calculation, the analysis we have made in this paper will serve as a basis for further study of the 
Hofstadter model through the elliptic curve techniques in algebraic geometry. This approach would 
allow us to follow a similar path as in the rational spectral curve case for the study of the spectrum 
problem of the Hofstadter model. This rich structure requires further study, and a scheme along 
this line of interpretation is under current investigation. Indeed, we hope that our efforts would 
eventually shed new light on the role of algebraic geometry in exactly solvable integrable models. 
In this paper we restrict our attention only to certain Hofstadter-like models, and leave possible 
generalizations and applications to future work. 
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